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We propose a new type of the Yang-Baxter equation (YBE) and the decorated Yang- 
Baxter equation (DYBE). Those relations for the fermionic i?-operator were introduced 
recently as a tool to treat the integrability of the fermion models. Using the YBE and 
the DYBE for the XX fermion model, we construct the fermionic i?-operator for the one- 
dimensional (ID) Hubbard model. It gives another proof of the integrability of the ID 
Hubbard model. Furthermore a new approach to the 5*0 (4) symmetry of the ID Hubbard 
model is discussed. 
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1 Introduction 

There have been much interests in strongly correlated electron systems. Among the exactly 
solvable ones, the one- dimensional (ID) Hubbard model, 

N ^11 

^ = - E E + sl^i+i-) + ^ EKt - ^)i^n - 2)- (1-1) 
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is the most interesting model. Here Cj^ and Cj^ are the fermion creation and annihilation 
operators and rija is the fermion number operator, 

%^ = sVj-- (1-2) 

The parameter U is the coupling constant describing the Coulomb interaction between 
electrons with opposite spins on the same site. Licb and Wu [1] diagonalized the Hamilto- 
nian (1.1) by means of the coordinate Bethe ansatz method under the periodic boundary 
condition (PBC), 

CN+l,a = Cjv+l,a = Ci^, (7 =ti ■ (1-3) 

The bulk properties of the ID Hubbard model have been investigated by analyzing the 
associated Bethe ansatz equations. [2] 

The integrability of the ID Hubbard model has been studied by Shastry [3-5] and 
Wadati et ai. [6-8] Shastry introduced the Jordan- Wigner transformation, [9] 

i-i 

CjT =exp|i7r^(nfc| - 

k=l 

i-1 N 
Cji =exp|i7r J^(nik; - 1)| expjiTT ^(nfe| - 1)}tj", (1.4) 

k=l k=l 

to change the fermionic Hamiltonian (1.1) into an equivalent coupled spin model, 

N N U ^ 

H = Y^^^ti^J + ^t^J^^ + + ^Z^^+i) + 1 E ^l^h (1-5) 
j=i j=i j=i 

where a and r are two species of the Pauli matrices commuting each other, and 

^f = ^K±i^J)' = ^(^;±i^/)- (1-6) 
Shastry constructed the L-operator and i?-matrix which satisfy the Yang-Baxter relation, 

R\2{ui,U2)Lal{ui)La2{u2) = -^a2(li2)-C'al (mi)-Ri2 (^1, ^2) ■ (1-7) 

In ref . [5] , the decorated Yang-Baxter equation (DYBE) for the free-fermion model was used 
to prove the Yang-Baxter relation for the coupled spin model, (1.7). It is a relation similar 
to, but not equivalent to the Yang-Baxter equation (YBE), Recently Shiroishi and Wadati 
[10, 11] found that the YBE and the DYBE for the free-fermion model can be generalized to 
a larger set of relations among the L-operators called tetrahedral Zamolodchikov algebra 
(TZA). [12] By use of the TZA, the Yang-Baxter equation for the i?-matrix, 

-Ri2(mi, M2)-Ri3(«1, U3)R23{u2, U3) = -R23(«2, M3)-Rl3(«l, UsjE^iui, U2) (1.8) 
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can be proved. [10, 11, 13] 

The coupled spin model (1.5) is often referred to as the ID Hubbard model, since it 
is related to the fermionic Hamiltonian (1.1) through the Jordan-Wigner transformation 
(1.4). However, because of the non-locality of the Jordan-Wigner transformation (1.4), 
the PBC for the fermion model (1.3) does not correspond to the PBC for the coupled 
spin model. Moreover the 5'0(4) symmetry [14-18], which is one of the most important 
properties of the Hamiltonian (1.1), is not clear in the coupled spin model (1.5). 

To discuss the fermionic Hamiltonian directly, we can employ the fermionic formulation 
of the Yang-Baxter relation (graded Yang-Baxter relation). [7] In this formulation, the 
transfer matrix that respects the PBC (1.3) can be constructed by taking the supertrace 
of the monodromy matrix. Recently it has been shown that the transfer matrix naturally 
enjoys the 5*0(4) symmetry. [19,20] Based on this fermionic formulation, the algebraic 
Bethe ansatz method has been developed by Ramos and Martins. [21,22] 

The aim of this paper is to present yet another approach to discuss the integrability 
of the ID Hubbard model using the fermionic ii!-operator. The notion of the fermionic 
i?-operator was introduced in relation to the XXZ fermion model. [23, 24] The fermionic 
i?-operator consists of the fermion operators and is a solution of the YBE. 

In this paper, we first consider the XY Z fermion model. The fermionic /2-operator 
for the XY Z fermion model, which satisfies the YBE, is explicitly obtained. The YBE 
ensures the commutativity of the transfer operators. The logarithmic derivative of the 
transfer operator yields the Hamiltonian of the XYZ fermion model under the PBC. 

Next we consider the DYBE for the fermionic i?-operator and find that the fermionic 
i?-operator related to the XY fermion model fulfills the DYBE. Combining the YBE and 
the DYBE of the fermionic i?-operator for the XX fermion models, we construct a new 
fermionic it!-operator and an L-operator for the ID Hubbard model. The transfer operator 
constructed from the L-operator generates the Hamiltonian and other conserved operators 
under the PBC. 

We also prove the YBE for the new integrable i?-operator by considering the fermionic 
variant of the tetrahedral Zamolodchikov algebra (TZA). Then we consider the generalized 
transfer operator which is constructed from the fermionic i?-operator. A new fermionic 
Hamiltonian is obtained from the transfer operator. 

It is possible to discuss the SO {4) symmetry of the ID Hubbard model using the 
fermionic -R-operator. The (anti) commutation relations between the fermionic i?-operator 
and the generators of the 5*0(4) algebra are obtained. Those relations are extended to the 
monodromy operator and the transfer operator. We can show the 5*0(4) invariance of the 
local conserved operators, though the transfer operator itself is not 5*0(4) invariant. 

The outline of this paper is as follows. In §2, we study the Yang-Baxter equation (YBE) 
for the XY Z fermion model. Fermionic /^-operator related to the XY Z fermion model is 
explicitly shown by use of Baxter's parameterization. In §3, the decorated Yang-Baxter 
equation (DYBE) for the fermionic /^-operator is introduced. It is necessary to impose the 
free-fermion condition on the Boltzmann weights in order that the DYBE is fulfilled. Some 
connections to the XY fermion model with an external field are also discussed. In §4, we 
discuss the integrability of the ID Hubbard model. The fermionic i?-operator for the ID 
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Hubbard model is constructed from the i?-operators related to the XX fermion models 
with spin up and down. The YBE and the DYBE arc utilized to show the Yang-Baxter 
relation, where a new L-operator is also introduced. In §5, we prove the YBE for the 
fermionic i?-operator of the Hubbard model. We use the fermionic variant of the TZA. 
The discussions proceed just in parallel with the case of the coupled spin model (1.5). We 
also use a generalized transfer operator which consists of the fermionic i?-operators. In §6, 
we give a new approach to the 5*0 (4) symmetry of the ID Hubbard model. The fermionic 
/?-opcrator has a natural 5*0(4) symmetry. We show that all the local fermionic conserved 
operators derived from the transfer operator are 5*0(4) invariant. The last section is 
devoted to the concluding remarks. 

2 Fermionic i?-Operator for the XY Z Fermion Model 

It is well understood that the Yang-Baxter equation (YBE) 

ni2{u - v)ni3{u)n23iv) = n2M'^Mni2{u - v), (2.1) 

plays the most important role in the quantum integrable models. [25-29] 
We look for a solution of the YBE (2.1) in the form 

Tljk{u) =a{u){-njnk + (1 - nj){l - n^)} + b{u){nj{l - n^) + (1 - nj)nk} 

+c{u){c]ck + clcj) - d{u){c]cl + CjCk), (2.2) 

where and Cj are the fermion creation annihilation operators satisfying the canonical 
anti-commutation relations 

{cj, Ck} = {cJ, 4} = 0, {c], Ck} = 5jk, (2.3) 
and Uj = CjCj. 

Substituting (2.2) into (2.1) and using (2.3), we get the following relations among the 
functions a (it), b{u), c{u) and d{u), 

a{u — v)c{u)a{v) + d{u — v)a{u)d{v) = h{u — v)c{u)b{v) + c{u — v)a{u)c{v), 
a{u — v)b{u)c{v) + d{u — v)d{u)b{v) = b{u — v)a{u)c{v) + c{u — v)c{u)b{v), 
c{u — v)b{u)a{v) + b{u — v)d{u)d{v) = c{u — v)a{u)b{v) + 6(m — v)c{u)c{v), 
a{u — v)d(u)b(v) + d{u — v)b{u)c{v) = b{u — v)d{u)a{v) + c{u — v)b{u)d{v), 
a{u — v)a{u)d{y) + d{u — v)c{u)a{v) = b{u — v)b{u)d{y) + c{u — v)d{u)a{v), 
d{u — v)a{u)a{v) + a{u — v)c{u)d{v) = d{u — v)b{u)b{v) + a{u — v)d{u)c{v). (2.4) 

These are nothing but the relations among the Boltzmann weights of the 8- vertex model 
(Baxter model). [25] Then, using Baxter's parameterization, we can solve (2.4) as follows, 

a{u) = sn{u + 2ri), b{u) = sn(-u), c(-u) = sn(2?7), 

d{u) = ka{u)b{u)c{u) = ksn{2rj) sii{u) sn{u + 2f]), (2.5) 
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where snu — sn{u, k) is Jacobi's elliptic function with the modulus k and r) is the anisotropy 
parameter. 

Thus we have found a new type of solution for the YBE (2.1). We call (2.2) the fermionic 
ii-operator (related to the XYZ fermion model). We remark that Destri and Segalini [23] 
have discussed the special case {k — 0) of (2.2). The fermionic i?-operator (2.2) can also 
be obtained from the i?-matrix of the XYZ spin model (Appendix A) . 

It is remarkable that TZjk{u) has a property, 

n,k{0) = sn2r]V,k, (2.6) 
where Vjk is the permutation operator, 

Vjk 1 - (c] - 4)(cj - Cfc) = 1 - Uj - rife + c]ck + cJ.Cj, 

This property is useful when we derive a related Hamiltonian from the transfer operator. 

Unlike the i?-matrix of the usual XY Z spin model, the fermionic i?-operator is not 
symmetric, 7^i2(m) 7^ 7^2i(w). Instead, the following relation holds 

n^2{u]k)^n2i{u]-k), (2.8) 

where we write the modulus k dependence explicitly. The fermionic i?-operator (2.2) also 
has the unitarity relation in the form 

7^l2(M)7^2l(-'u) = {-SYi^u + SYi^2ri)l, (2.9) 

as well as the property 

[7^,•fc(M), (1 - 2n,)(l - 2nfc)] = 0. (2.10) 

Now we introduce the monodromy operator 

%{u)^naN{u)...naM. (2-11) 

where the suffix a means the auxiliary space. The transfer operator is given by the "super- 
trace" of the monodromy operator 

T{u)^^iYaTa{u). (2.12) 

Here we define the supertrace Str as 

Str„X = „(0|X|0)„ - a{l\X\l)a, (2.13) 
with the auxiliary fermion Fock space being 

c„|0)„ = 0, |1)„ = 4|0)„, „(0| = (|0)„)t, a(l| = a(0|c„, (2.14) 
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„(0|0)„ = „(1|1)„ = 1. (2.15) 
When we write the monodromy operator Ta{u) as 

%{u) = A{u){l - Ua) + B{u)ca + C{u)c\ + L>(ii)n„, (2.16) 
the transfer operator (2.12) is given by 

t{u)^ A{u)- D{u). (2.17) 

The YBE for the fermionic i?-operator leads to the global Yang-Baxter relation for the 
monodromy operator (2.11), 

ni2{u - v)Ti{u)T2{v) = T2{v)Mu)nu{u - v). (2.18) 

Rewriting (2.18) as 

TZuiu - v)T,{u)T2{v)n^^\u -v)= T2{v)T,{u), (2.19) 

and taking the supertrace (2.13) of both sides, we can show the commutativity of the 
transfer operators 

[riu),riv)]^0, (2.20) 

which ensures the existence of enough number of the fermionic conserved operators. In 
this way we can establish the integrability of the fermion model by means of the fermionic 
i?-operator. 

Using (2.6) in (2.12) with (2.11), we have 

r(0) = {sn2r]f Str,(P„^ . . . Pa2Vai) 

= (sn 27;)^ • • • rN-i,N: (2.21) 

where we have used some fundamental properties of the permutation operators 

Str„P„, = 1. (2.22) 
Thus we find that t(0) is proportional to the (left) shift operator U for fermions, 

U = Vi2V2z...rN-i,N. (2.23) 
It is easy to confirm that the shift operator U acts on the fermion operators as 

Ucj = Cj+iU, j = 1, . . . , - 1, 

Ucn = ciU. (2.24) 



In a similar way as the spin models, the Hamiltonian of the XYZ fermion model is 
derived from the logarithmic derivative of the transfer operator 



du 

N 



u=Q 



E«i;i?.. (2-25) 



where the two-point Hamiltonian density is given by 

+ 6'(0)(c]c,+i + c]+i9) + d'mc]c]^, - 99+1)}. (2.26) 

Here and hereafter, ' means the derivative with respect to the spectral parameter 

a'(0) = cn277 dn27;, 6'(0) = 1, d'{Q) = k sn^2r}. (2.27) 

In (2.25), we have used a fact that t(0)~^ is proportional to the shift operator. In terms 
of an equivalent i?-operator TZjk{u) = Vjk'R-jkiu) (c.f. Appendix A), the two-point Hamil- 
tonian density is expressed as 



(2.28) 

u=0 



One may wonder if our discussion using the fermionic i?-operator goes fully in parallel 
with the standard treatment of the integrable spin models and no new aspect appears. In 
fact we can proceed the discussion almost in the same way as the spin models. However 
we like to point out an important difference between the spin model and the fermion 
model here. While the Hamiltonian (2.26) is hermitian, it is not invariant under the space 
inversion 

j^N-j + 1. (2.29) 

It is easy to see that the modulus k in the Hamiltonian (2.26) changes its sign under the 
transformation (2.29). This point makes a difference between the XYZ spin model and 
the XYZ fermion model. Actually it is a consequence of the property (2.8). 
In the case of A; = 0, the fermionic i?-operator (2.2) reduces to the one with 

a{u) = sm{u + 2r)) h{u) = sinu c{u) = sm2r] d{u) = 0, (2.30) 

which we have used in ref. [24] to discuss the integrability of the XXZ fermion model 

1 ^ 

n''''^ =— ^ Yl {S^^i+i + sViC.- + cos 277 (n,n,+i + (1 - n,)(l - n,+i)) }. (2.31) 
' i=l 
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In this case, the fermionic i?-operator becomes symmetric and the Hamiltonian (2.31) is 
invariant under the space inversion (2.29). 

We note that the periodic boundary condition (PBC) for the fermion operators is 
satisfied in (2.26), 

cn+1 = ci, cjy^i = c\. (2.32) 

It is sometimes necessary to consider the anti-periodic boundary condition (APBC) for the 
fermion operators 

cn+1 = -ci, cjv_^i = -c\. (2.33) 
To explain this, we consider the following twisted monodromy operator 

fa{u) = (1 - 2na)Ta{u) 

^{l-2na)naN{u)...nal{u). (2.34) 

Then, thanks to the property (2.10), the global Yang-Baxter relation for the twisted mon- 
odromy operator holds, 

7^l2(n - v)f,{u)f2iv) = f2iv)f,{u)nuiu - v), (2.35) 

which leads to the commutativity of the twisted transfer operator f{u) — StYaXiiu), 

[f{u),f{v)]^0. (2.36) 

We can show that the twisted transfer operator generates the Hamiltonian (2.26) under 
the APBC (2.33). In fact the logarithmic derivative of f{u) yields a Hamiltonian 



au 



u=0 



N-1 ^ 

- 6'(0)(4ci + c\cn) - c?'(0)(44 - cjvci)}. (2.37) 

Here we have used a relation 

(1 - 2n^) Str, {(1 - 2na)VaNn^/''} = (1 - 2n;v)^m ^(1 - 2njv), (2.38) 
as well as the identities 

(1 - 2niv)civ(l - 2n^) = -cn, (1 - 2nAr)4(l - 2n^) = -4. (2.39) 

Note that if we express Ta{u) as (2.16), the twisted transfer operator takes the following 
form, [23] 

f{u)=A{u) + D{u). (2.40) 



3 Decorated Yang-Baxter Equation 

The decorated Yang-Baxter equation (DYBE) is an algebraic relation similar to the YBE 
(2.1), but not equivalent to the YBE. Shastry [5] showed that the free-fermion condition 
on the Boltzmann weights is necessary so that the solution of the YBE also satisfies the 
DYBE. [5,11] The DYBE is considered to be a hidden algebraic structure of the free-fermion 
model (equivalently, XY model). 

We discuss the DYBE in the context of the fermionic i^-operator. The DYBE for the 
fermionic i?-operator is defined by 

TZuiu + v){2ni - l)nn{u)n23{v) = 7^23(^;)7^l3(^i)(2nl - l)7^l2(^^ + v). (3.1) 

The DYBE (3.1) looks alike the YBE (2.1), but the operator (2ni - 1) inserted in both 
sides can not be absorbed by the redefinition of Rjk{u). It is possible to simphfy the DYBE 

by use of the YBE as follows. 

First, by multiplying (2ni — l)(2n2 — 1) to the YBE (2.1) from the left and (2^3 — 1) 
from the right, we get a relation 

7^l2(^^ - v){2ni - I)ni3{u){2n2 - l)7^23(v)(2r^3 - 1) 

= (2712 - l)7^23(^^)(2r^3 - l)7^l3(^^)(2r^l - l)nu{u - v), (3.2) 
where we have used a general property (2.10) and an equivalent relation, 

(2n,- - l)njk{u){2nk - 1) = {2nk - l)7^,•fc(^^)(2r^,■ - 1). (3.3) 
Then comparing cq. (3.2) with the definition of the DYBE (3.1), we find a condition 

{2nj - l)njk{u){2nk - 1) = 7^,•fe(-^^). (3.4) 

Hence we see that the DYBE is equivalent to the condition (3.4) when we have the YBE. 
In terms of the Boltzmann weights, the condition (3.4) is rephrased as 

a{u) = a{—u), —b{u) = b{—u), c{u) = c{—u), —d{u) = d{—u). (3.5) 

Using the explicit parameterization (2.5), we find that the condition (3.5) corresponds to 

2ri = K, 3K, (3.6) 

where K is the complete elliptic integral of the first kind. In this case, the free-fermion 
condition 

a^iu) + b^u) = c\u) + (fiu), (3.7) 

is fulfilled. Thus as in the case of the spin models, we have shown that the free-fermion 
condition is required for the validity of the DYBE (3.1). 
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Hereafter we restrict our consideration to the case 2r) — K. Then, from (2.5), the 
Boltzmann weights a{u), b{u), c{u) and d{u) are exphcitly given as 

a{u) = cdu, h{u) = snu, c{u) = 1, 

d{u) = kcdu snu, (3.8) 
where cdu = cnu/ dnu. The corresponding Hamiltonian is the XY fermion model 

N 

Ti^Y ^ J2 {c]cJ+^ + c]+,9 + k (c]c]+, - c,c,+i) } . (3.9) 
i=i 

In the special case A; = 0, this Hamiltonian reduces to the XX fermion model 

N 

^'''' = E(sVi + 4+i9)- (3-10) 

i=i 

As an interesting application of the YBE (2.1) and the DYBE (3.1), we shall prove the 
integrability of the XY fermion model in an external field, 

N / 1\ 

^xym = Yl {4^J+^ + + k (c]c]+i - c,c,+i) +i,(nj--]]. (3.11) 

j=i V / 

The external field fi can be regarded as the chemical potential for the XY fermion model. 
Since the term ^^^(nj — 1/2) does not commute with the Hamiltonian (3.9) unless A; = 0, 
the addition of the extra term is not trivial at all. It is well known that the model 
(3.11) can be diagonalized by means of the Fourier transformation through the Bogoliubov 
transformation. [30] The case k = ±1 corresponds to the transverse Ising model. [31] 

The following argument is in common with the one explored for the spin models. [3,11] 
Prom now on, let 7ljk{u) denote the fermionic i?-operator for the XY fermion model. 
Taking a linear combination of the YBE (2.1) and the DYBE (3.1), we have 

^ani2{u-v)+(5ni2{u + v){2ni - l)^ni^{u)n2z{v) 

= 7^23(^;)7^l3(^i){Q;7^l2(^^ - v) + (5{2ni - l)7^l2(^^ + ^;)}, (3.12) 

where a and (3 are not operators but constants. 

Now we look for a solution of the Yang-Baxter relation 

nl2{u,v)Ci^{u)C2^{v) ^ C2^{v)Ci^{u)nl2{u,v), (3.13) 

assuming the i?-operator TZl2{u,v) and the L-operator Cjk{u) in the form, 

n'^^(u,v) = a7^l2 {u-v)+ /57^l2 (u + v) (2ni - 1) , (3. 14) 

Cjk{u) = TZjkiu) exp{h{u){2nj - 1)} . (3.15) 
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Here a, (5 and h{u) are functions of the spectral parameters to be specified later. 
Comparing eq. (3.13) with eq. (3.12), we get a relation 

Xi{u)X2{v)nl2{u,v)Xi{u)-^X2{vy^ = ani2{u - v) + (5{2ni - l)7^l2(^^ + 1;), (3.16) 
where 

Xj{u) = exp{h{u){2nj - 1)}. (3.17) 

For k ^ (that is, the genuine XY fermion model), the relation (3.16) gives two 
conditions, 

f3 



a 



teinh{h{u) - h{v)) (3.18) 



^'^1'' ^ '"I = teinh(h(v) + h(u)). (3.19) 
ad{u — V) 

Equations (3.18) and (3.19) give the ratio of a to /3 and constraints on h{u) and h{v). 
Substituting (3.8) into (3.19), we find that the constraints are given by 

smh.2h(u) smh.2h(v) 

^ = ^ = constant, 3.20 

sc 2u sc 2v 

where scu — snu/ cnu. In conclusion, we have proved that the fermionic i?-operator 

n'^^{u,v) = 7^12(M -w) +tanh(/i(M) - h{v))nuiu + v){2ni - 1), (3.21) 

satisfies the Yang-Baxter relation (3.13) with the L-operators (3.15) under the constraints 
(3.20) . 

Prom eq. (3.13) we get the global Yang-Baxter relation for the monodromy operator 
Xa{u)=LaN{u)...L,M. (3.22) 

as 

TV[^{u,v)r^{u)%{v) = r2{v)r^{u)'R}i^{u,v), (3.23) 

which leads to a one-parameter family of the commuting transfer operators, 

[r(«),r(t;)] = 0, 

r(M) = Str„T„(M). (3.24) 

If we identify the constant in eq. (3.20) with /u/2, we recover the Hamiltonian (3.11) 
by the logarithmic derivative of the transfer operator 

au Im=o 

+ (3.25) 

i=i ^ ^ 

Here we have used the relations /i(0) = 0, /i'(0) = ///2. 
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4 Fermionic i^-operator for the ID Hubbard model 

The Hamiltonian of the ID Hubbard model (1.1) consists of two XX fermion models (3.10) 
for up and down spins with an interaction term between them. The fermionic i?-operator 
without an interaction term is given by 

njk{u) = n'^^{u)n^i{u), (4.i) 

where TZ^'^J {a =ti) denote the fermionic i?-operators for the XX fermion model, i.e., 

a{u) — cosu, b{u) — sin 14, c{u) — 1, d{u) — 0. (4.2) 

Since both 7^^|.^(m) and Tl^j\^{u) satisfy the YBE and the DYBE, the product 'R.jk{u) also 
satisfies the YBE 

ni2{u - v)ni3{u)n23{v) = n23{v)ni3{u)ni2{u - v) (4.3) 

and the DYBE 

'JZi2{u + v){2ni^ - l)(2ni| - l)7ei3(M)^23(t') 

= n23{v)ni3{u){2n^ - l)(2ni; - l)nu{u + v). (4.4) 

A hnear combination of (4.3) and (4.4) yields 

^^aUuiu -v) + pTZuiu + v){2ni^ - l){2nii - I)}ni3{u)n23{v) 

[a1Zi2iu -v)+ /3(2niT - l)(2nu - l)1Z^2iu + i;)}. (4.5) 

For a moment, a and f3 are arbitrary. As in the case of the XY fermion model in an 
external field, we look for a solution of the Yang-Baxter relation, 

n\^{u,v)Ci3{u)C23{v) = C23{v)C,3{u)n\^{u,v) (4.6) 

in the form 

T^u{u,v) = a1Zi2{u-v)+ p1Zi2{u + v){2ni^ - l){2n^ - 1), (4.7) 
Cjk{u) ^Hjkiu) exp {-h{u){2nj^ - l)(2n,-^ - 1)}. (4.8) 

Comparing eq. (4.5) with the Yang-Baxter relation (4.6), we get a relation 



Ti {u)l2 {v)TZ\^ {u, v)Ii (u) -^Mv)-^ 

= aUuiu -v)+ /3(2n^ - l)(2ni^ - l)1Zi2{u + v), (4.9) 
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where 

X,{u) = exp{-/i(«)(2n,T - l)(2n,^ - 1)}. (4.10) 
Prom (4.9), we obtain two conditions 

—7 { = - tanh(h(u) - h(v)), f = - tanh(h(u) + hiv)). (4.11) 

aa{u — v) ah{u — v) 

which give the ratio of o; to /3 and constraints on h{u) and h{y). The constraints are given 
more exphcitly by 

sinh2/i(-u) sinh2/i(f) U /a ^r,\ 

sin 2u sin 2v 4 

To sum up, we have obtained the fermionic it!-operator for the ID Hubbard model 

cos[u + v) 

X 7^^5(^i + v)n^i{u + v){2n^ - l)(2ni| - 1), (4.13) 

which satisfies the Yang-Baxter relation (4.6). One of the remarkable properties of the 
fermionic i?-operator (4.13) is that its dependence on the the spectral parameters is not a 
difference type u — v. We call it the "non-additive" property. The non-additive property of 
the fermionic i?-operator (4.13) allows us to generalize the Hamiltonian of the ID Hubbard 
model (see §5). 

From (4.6) we get the global Yang-Baxter relation, 

n\^{u,v)T^{u)%{v) ^%{v)T,{u)n\^{u,v) (4.14) 

for the monodromy operator 

Ta{u)^CaN{u)...Cal{u). (4.15) 

Now we define the supertrace Str** for the ID Hubbard model as 

Str^ X = „(0|X|0)„ - „(t \X\ t)„ - ail \X\ Da + a(it |^| ti)a, (4-16) 

where the auxiliary fermion Fock space is given by 

Ca.|0)a = 0, I T)a = 4|0)., | i)a = cl\0)a, \ TDa = cl^cljO)„, 
a(0| = (|0)„)t, „(t I = „(0|C„T, a{i\ = a{0\Cai, a{iU = a{0\CalCaV (4-17) 

Then the transfer operator of the ID Hubbard model is defined by 

T{u) = Str'^Uu). (4.18) 
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Prom the global Yang-Baxter relation (4.14), we can show the existence of a commuting 
family of the transfer operators 



[t{u),t{v)] = 0, 



(4.19) 



which proves the integrability of the model. 

Using the relations h{0) — and h'{0) = U /A, we can obtain the Hamiltonian of the ID 
Hubbard model as a logarithmic derivative of the transfer operator, 



where the PBC for the fermion operators (1-3) is satisfied. 

In this way we have proved the integrability of the ID Hubbard model using the 
fermionic it!-operator. It should be emphasized that we have not used the Jordan- Wigner 
transformation at all, which changes the boundary condition and the symmetry of the 
model. [20] We shall see in §6 that the fermionic i?-operator (4.13) naturally enjoys the 
5*0(4) symmetry. 

5 Yang-Baxter Equation for the New Fermionic R- 
Operator 

In the previous section we have shown that the fermionic i?-operator (4.13) satisfies the 
Yang-Baxter relation (4.6) with the L-operator (4.8). It is natural to expect that the 
fermionic i?-operator itself fulfills the Yang-Baxter equation (YBE), 



In the case of the coupled spin model, the YBE of the i?-matrix was proved in ref. [10] 
using the tetrahedral Zamolodchikov algebra (TZA). [12,13] In this section, we introduce 
the TZA for the fermionic i?-operator and proves the YBE (5.1). 

The TZA is defined by the following set of relations among operators C^^^ and C}j^, 

^12^\z^2d, = 5^ Sde}^2Z^\z^l2^ O, . . . , / = 0, 1, (5.2) 




= -EE^S? + TE(2-.T-l)(2n,,-l), 



(4.20) 



j=i <T=n 3=1 



n\^{Ui, U2)n\^{Ui, M3)7^23(M2, «3) = ^23(^2, M3)7^^3(Ml, U^)n\^{Ui, U2). 



(5.1) 



d,e,/=0,l 




(5.3) 
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where 'R}'^^{u) and TZ^^^{u) are the fermionic it!-operators for the XX fermion model as 
before. Then we have found the following relations which give the TZA (5.2), 

^O(a) ^O(a) ^O(a) _ ^O(a) ^O(a) ^O(a) ^0{a) ^l{a) ^l{a) _ ^l(a) ^l(a) ^0{a) ,^ 

''-12 ''-IS ''-23 ~ ''-23 ''-13 ''-12 ) ''-12 ''-13 ''-23 ~ ''-23 ''-13 ''-12 W-^J 
^l(a)^l{a)^0(a) _ Ma) ^1(<7) ^0(<7) ^1(<7) _ ^l(a) ^O(a) ^1(<t) ,^ 

'-12 ''-13 ''-23 ~ ''-23 ''-13 ''-12 ; ''-12 ''-13 ''-23 ~ ''-23 ''-13 ''-12 ? \^-^) 
pl{a) pl(a) _ qlll pO{a) pO{u) qUI ^O(cr) ^l(o-) ^O(cr) , qIII a.0(o-) ^O(o-) ^1((t) 

''-12 ''-13 ''-23 — *-'001''-23 ''-13 ''-12 "r *-'010''-23 ''-13 ''-12 "r *-'l00''-23 ''-13 ''-12 ) 

(5.6) 

pO{a) pO{a) pl{a) _ eOOl r^^r) pl{a) pl{a) , qOOI /-O(cr) ^O(o-) ^l(cr) , qOOI /-O(o-) ^1((t) ^0((t) 
''-12 ''-13 ''-23 ~ *-'lll''-23 ''-13 ''-12 *-'l00''-23 ''-13 ''-12 '-'010''-23 ''-13 ''-12 ' 

(5.7) 

^0(<t) ^1(<t) ^O(o-) _ cOlO /-IW /-IW /.l(cr) I cOlO /'0((t) ^O(ct) ^l(cr) qOIO /-IW /.O(cr) ^0(<t) 
''-12 ''-13 ''-23 ~ *-'lll''-23 ''-13 ''-12 "r *-'l00''-23 ''-13 ''-12 "r *-'001''-23 ''-13 ''-12 ) 



(5.8) 

) 

12 -^13 -^23 ~ '^111-^23 -^13 -^12 ^ "010-^23 -^13 -^12 ^ "001-^23 -^13 -^12 ' 



^l(o-) ^O(o-) ^0((t) _ clOO /-^o-) ^l(o-) ^l(o-) clOO /-OM ^l(cr) ^O(cr) clOO /-IM ^O(o-) ^O(cr) 
■'-IS ■'-93 — *-'l11^'-9S ■'-IS ■'-19 ~r )Jr|lfT'-9S ■'-IS ■'-19 "T Orini ■'-9S ■'-IS ■'- 



(5.9) 



where u =t or J,, and the coefficients S'^^ are given by 



_ sin(Mi + U2) cos(mi + M3) _ sin(iii + U2) sin(M2 + M3) 

'-'001 ~ „ /„, „. „. '-'010 



cos(iii — U2) sin(Mi — M3) ' cos(ui — U2) cos(m2 — u^) ' 

^111 _ sin(^2 + ^3) cos(^i + ^a) ^^^-^ _ sin(^i - ^2) cos(^i - ^3) 

cos(w2 — 1*3) sin(ui — U3) ' cos(wi + U2) sin(tii + u^) ' 

^001 _ sin(Mi - U2) sin(M2 + %) ^001 ^ sin(M2 + %) cos(mi - %) 

cos(hi + M2) cos(m2 — M3) ' cos(ii2 — u^) sin(Mi + -us) ' 

^010 _ sin(Mi - U2) sin(M2 - ^3) ^010 _ sin(Mi - M2) cos(mi + M3) 

^"^"^ COs(mi + U2) 008(^2 + M3) ' cos('Ui + U2) sin(-UiM3) ' 

^010 ^ sin(M2 - M3) cos(mi + M3) ^^qq ^ sin(M2 - M3) cos(mi - u^) 

cos(u2 + 1x3) sin(ui — ^3) ' cos(w2 + M3) sin(-Ui + M3) ' 

^^00 _ sin(tti + U2) cos(mi - u-i) ^^qq _ sin(-»i + -^2) sin(M2 - u^) , . 

'-'nin — / \ • / \ ' '-'001 — / \ / \ ■ 10. iU) 

cos(mi — U2) sin(iii + 1*3) cos(mi — U2) C0S\U2 + Uz) 

Note that eq. (5.4) and eq. (5.5) are equivalent to the YBE (4.3) and the DYBE (4.4) 
respectively. In this sense, the TZA (5.2) can be regarded as a generahzation of the YBE 
and the DYBE. 

Wc also remark that the products C^^^ C''^^^ £^^^ (a, b,c — 0, 1) are not linearly inde- 
pendent and satisfy the following relations, 

^0{a) ^O(a) ^0(a) _ ^0(a) ^l(a) ^l(a) ^,^1M ^0{a) ^l(a) ^l{a) ^l(a) ^0(a) (r^TU 

-'-12 '-IS -'-23 "•^''-12 -'-13 -'-23 i'-'-12 -'-13 -'-23 "'"^-'-12 -'-13 -'-23 ' ^--'--'-J 

^l{a)^l{a)^l(<7) _ /^1(<7)^0((7)^0(<7) , ^0{a) ^l{a) ^O(a) , ^0{a) ^0{a) ^l{a) 

'-12 ''-13 ''-23 "•^■'-12 ''-13 ''-23 T" i/ ■'-12 ''-13 -'-23 "'"^-'-12 ''-13 ''-23 ' W'-'-^J 
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where 



sin(ui — Us) cos(w2 — Us) cos(iii — U2) cos(u2 — 1*3) 



sin(ui + Us) cos(ii2 + Us) ' cos(iii + U2) cos(w2 + Us) ' 

cos(iii — 1^2) cos(iii — Us) I sin(wi + us) cos (1*2 + us) 
cos(wi + M2) cos(iti + Us) ' sin(-ui — -^3) cos(-U2 — ""3) ' 

, _ cos(mi + U2) cos(m2 + us) , _ sin(Mi + us) cos(mi + M2) 1 Q^ 

^ cos{ui — U2) cos{u2 — UsY sm{ui — Us) cos{ui — U2) 

These relations (5.4) - (5.13) have been verified using mathematica. 

Now taking into account the form of the fermionic i?-operator (4.13), we look for a 
solution of the YBE (5.1) in the form 

T^)kiuj, Uk) = Cf^\uj, Uk)Cfi\uj, Uk) + ajkCf^\uj, Uk)Cff}\uj, Uk). (5.14) 

We substitute eq. (5.14) into the YBE (5.1). Using the TZA (5.2) and the relations (5.11) 
and (5.12), we obtain the conditions 

ais sin 2(mi + us) = sin 2(-Ui + U2) + a2s sin 2(^2 + Us) 

— sin 2{ui — U2) sin 2{ui — us)- (5.15) 

ai2 0(is 

If we assume aij in the form 

— — ta,nh(h(ui) — hiuh)). (5.16) 

COs{Uj + Uk) 



and impose the constraints 

J = 1,2,3, (5.17) 



sinh2//((/j) _ U 



sin 2uj 4 ' 

then one can easily find that the conditions (5.15) are satisfied. In conclusion, the fermionic 
i?-operator for the ID Hubbard model, 

T^%iuj,Uk) =^5?(«,- - Uk)n'^^{uj - Uk) - ^^^^!^tanh{h{uj) - h{uk)) 

X n^^^ {uj + Uk)nf,^ {u, + Uk) (2n,^ - 1) (2n,-^ - 1) , (5. 18) 

is a solution of the YBE (5.1) under the constraints (5.17). 

Besides the YBE (5.1), the fermionic /^-operator (5.18) has the following properties 

^'-''^ 

n%{uo,Uo)=V,k, (5.20) 
Uk)nlj{uk, Uj) = p{uj, Uk)l, (5.21) 
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where 

p{uj, Uk) = cos^{uj — Uk) {cos^{uj — Uk) — taji\i'{h{uj) — h{uk)) cos^{uj + u^)^ . 

(5.22) 

Here the permutation operator is defined by 

T? _ 'r)(T)-p(i) 

Vf^ = 1 - {cl - cl){c,^ - Ck.), a =n . (5.23) 

Due to eq. (5.19), we can recover the Yang-Baxter relation (4.6) by putting = in the 
YBE (5.1). 

Using the /^-operator (5.18), we can introduce an inhomogeneous model as 

%{u, {uj}) = 7^^Jv(^^,^^iv) ■ ..n%{u,Uj) . . .n'^,{u,ui), (5.24) 
where uj {j = 1, . . . , N) are the inhomogeneous parameters obeying the constraints 
sinh2h{uj) U 



sin 2uj 4 ' 



j = l,...,N. (5.25) 



Then from the YBE (5.1), we have the global Yang-Baxter relation 

n\^{u, v)T,{u, {uj})T2{v, {uj}) =%{v, {u^])r^{u, {u^])n\^{u, v), (5.26) 
which leads to the commutativity 

[T(«,K}),r(^;,K})] = 0. (5.27) 
Because of the non-additive property of the i?-operator (5.18), the simplest choice 

Ta{u,uo) = 7^^^(^^,^^o) . ..Tl^j{u,uo) . ..Tl^^{u,uo) (5.28) 

still generalizes cq. (4.15). Actually we can derive a fermionic integrable Hamiltonian which 
generalize the ID Hubbard model. Note that t{uq,uo) is the shift operator regardless of 
the parameter uq, 

r{uo, uo) = U 

= Vi2V23...Vn-i,n. (5.29) 
Then the logarithmic derivative of the generalized transfer operator yields a new fermionic 
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Hamiltonian as follows, 



71 = - t{uo,Uo) 



-1^ 
du 



t{u,Uo) 



u=uo 



N 



+ 



U 



4cosh2M«;) ^ { ^°«'^o(2n,T - 1) - sin^ ^^o(2r^,+lT - 1) 
j=i 

+ sin 2ko (c]+^Cjt - cJ^Cj+it)| 

X I cos^ uo{2nji — 1) — sin^ Mo(2nj+ij — 1) 

+ sin 2mo(c]+i^Cj; - cj^^+i J } . (5.30) 

If we take = or 7r/2, the Hamiltonian (5.30) reduces to the ID Hubbard model (1.1). 
The local higher conserved operators can also be obtained by the formula 



d" 



du 



-ln(r(Mo,Mo) V(m,mo)) 



n=l,2,... , 



(5.31) 



U=Uo 



where Ii = —Ti. The meaning of the insertion of the term r(uo,Uo)~^ in (5.31) will be 
clarified when we discuss the 5*0(4) symmetry in §6. 

Remark that in terms of TZjkiu^uo) = Vjk'R-jkiUjUo) we can express the Hamiltonian 
(5.30) as 



N 



= - X^ ^TZjkiu, uo) 



U=Uo 



(5.32) 



The higher conserved operators are also expressed in terms of TZjk{u,Uo), though it is 
very difficult to get their explicit forms. 



6 SO (4) Symmetry 

The S0{4:){= [SU{2) x SU{2)]/Z2) symmetry is one of the most important properties of 
the Hubbard Hamiltonian (1.1) [14-18]. The two SU{2) come from the spin SU{2) 

N N ^ N 

j=i j=i j=i 

and the charge SU{2) (r^-pairing SU{2)) 

N N ^ N 

i=i j=i j=i 
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We assume that the lattice has an even number of sites under the PBC (1.3). Then the 
six generators above commute with the Hamiltonian (1.1). 

In ref. [19,20], the 5*0(4) symmetry of the fermionic transfer matrix was discussed based 
on the fermionic formulation by Olmedilla et ai. [7] In the following we discuss the 5*0(4) 
symmetry of the fermionic i?-operator and the generalized transfer operator (5.28) using 
the similar discussion (c.f. ref. [32]). For this purpose, we introduce the local generators 
of the 5*0(4) algebra as follows. 



Si 



^1 



(6.3) 



The SU{2) symmetry of the fermionic /^-operator for the ID Hubbard model (5.18) is 
represented by the following commutation relations. 



n%{u,v),s^ + s^ 
n%{u,v),s^ + si 



= 0, 
= 0. 



= 0, 



(6.4) 



It is easy to see that the relations (6.4) extend to the SU{2) symmetry of the monodromy 
operator 

Ta{u, Uo) = Tl^Ni^, uo)... Ti\^{u, Uo) 



as 



%{u,uo),s: + s^ 



0, 



Taking the supertrace Str^ of (6.5), we obtain 



0, 



(6.5) 



(6.6) 



which shows that the transfer operator t{u,Uq) is the spin SU{2) invariant. 

Next we consider the charge SU{2) symmetry. For the generator r/J, the same relation 
as (6.4) holds. 



TZ%{u,v),r)l + r)l 



-0, 



(6.7) 



which extends to the symmetry of the monodromy operator and the transfer operator as. 



0, 



0. 



(6. 



As for the other generators rjf, we do not have the relations like (6.4). Instead, we have 



discovered the following anti-commuting relation for rj^, 



(6.9) 
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The local relation (6.9) extends to the identities for the monodromy operator, 

r„K«o),r7^] + [Ta{u,uo),v^} = 0, (6.10) 
which can be proved as follows, 

= ^a7v(^> «o) • • • 7^a2(?^, Mo)^al(«, Uo){Va ' Vl + V2 ■ ■ ■ + Vn) 
= n^Niu, Uo)... Tla2i'^, Uo){-Va + vt + ■ ■ ■ + ^^)^al(^*> ""o) 



Here we have used the fact that N is even. 
Taking the supertrace of (6.10), we obtain 



[r{u,Uo),r]^^ = 0. 



(6.11) 



(6.12) 



Namely the transfer operator anti-commutes with the generators 77^. Thus contrary to the 
naive anticipation, the transfer operator itself does not commute with the generators of 
the charge SU{2) symmetry. However a combination t{uo,Uo)~^t{u,Uo) commutes with 



t{uo,uo) ^T{u,uo),ri^ 

due to the relation 

{r{uo,Uo),r)^^ = 0. 



= 0, 



(6.13) 



(6.14) 



Of course, t{uo, uq) ^t{u, uq) commutes with the other generators of the S'0(4) algebra. 
Since we have defined the local conserved operators by eq. (5.31), 



4 = ^In (r(Mo,Mo) V(m,mo)) 



n=l,2,... 



(6.15) 



u=uo 



it follows that they are all SO{A) invariant. In particular, we have proved that a gener- 
alized Hubbard model (5.30) has the 5*0(4) symmetry. It is also possible to prove the 
5*0(4) invariance of the Hamiltonian (5.30) directly using the expression (5.32) and the 
commutation relations 



n%{u,v),sf + s^ 



0, 



0, 



a — ±,z, 



0. 



(6.16) 
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where TZ^^{u,v) = VjkTl%{u,v) as before. 

The relation (6.14) means that the shift operator t) — t{uo,uo) anticummutes with r)^ 

Ur)^ = -r]^U. (6.17) 

If we define the lattice momentum operator ft [19, 32] by 

[/ = expin, (6.18) 
the relation (6.17) implies that 

tlr)^ - rj^tl = nr)^. (6.19) 

This indicates that the generators f]"^ changes the momentum eigenvalue by tt. [15,19] Note 
that the lattice momentum operator ft is defined in (6.18) mod 27r. 
The lattice momentum operator ft is also a conserved operator 



n,/n 



0, (6.20) 



but it is not contained in {In}- Actually ft is a non-local operator (we refer ref. [19] for 
the detailed discussions on ft). Furthermore the relation (6.19) suggests that it does not 
have the 5*0(4) symmetry. 



7 Concluding Remarks 

In this paper, we have presented a new fermionic approach for the integrability of the ID 
Hubbard model. The central object in our approach is the fermionic i?-operator, which 
we introduced recently to discuss the integrability of the XX Z fcrmion model. We have 
shown that the fermionic R- operator can be generalized to the XYZ fcrmion model. It 
gives an operator solution of the Yang-Baxter equation (YBE). We also considered the 
decorated Yang-Baxter equation (DYBE) for the fermionic i?-operator and found that the 
free-fermion condition is necessary. In other words, the fermionic i?-operator related to 
the XY fermion model fulfills the YBE and the DYBE. 

The ID Hubbard model can be regarded as two XX fermion models with an interaction 
term. Using the YBE and the DYBE for the XX fermion model, we have established a 
new kind of the Yang-Baxter relation for the ID Hubbard model. This gives an alternative 
proof for the integrability of the ID Hubbard model. We have obtained the fermionic R- 
operator for the ID Hubbard model and confirmed that it satisfies the YBE. The connection 
between the fermionic i?-operator and the graded Yang-Baxter relation for the ID Hubbard 
model [7] will be discussed in a separate paper. [33] 

From the fermionic /^-operator of the ID Hubbard model, we have constructed a new 
fermionic transfer operator. Because of the non-additive property of the i?-operator, the 
transfer operator generates a fermionic integrable Hamiltonian (5.30) which generalize the 
ID Hubbard model. 
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One of the advantages of our fermionic approach is that we can discuss the 5'0(4) 
symmetry at the level of the (local) fermionic J?-opcrator, which immediately leads to the 
5*0 (4) symmetry of all the local conserved operators derived from the transfer operator. It 
is known that the 5*0(4) symmetry of the ID Hubbard model is extended to the Yangian 
symmetry on the infinite lattice [34,35]. It is interesting to investigate the relation between 
the fermionic i?-operator and the Yangian symmetry. 

Acknowledgment 

The authors arc grateful to Y. Komori for useful discussions and careful reading of the 
manuscript. This work is in part supported by Grant-in- Aid for JSP J Fellows from the 
Ministry of Education, Science, Sports and Culture of Japan. 



APPENDIX A: A Simple Derivation of the Fermionic i?-Operator 

The it!-matrix of the XYZ spin models in terms of the Pauli spin matrix {a^', a~ and cr|) 
is given by 

Rjk{u) =a(M)ii^^ + h{uy-^^ + c{u){a^a^ + a^aj) 

+ d{u){ajat + ajal), (A.l) 

where the Boltzmann weights a{u), b{u), c{u) and d{u) are the same as (2.5). The i?-matrix 
(A.l) satisfies the Yang-Baxter equation 

Ri2iu - v)R,3{u)R23iv) = R23iv)RMR^2iu - v). (A.2) 

We like to obtain the fermionic i?-operator (2.2) from the i?-matrix (A.l) by means of the 
Jordan- Wigner transformation [30] 




(A.3) 



However it is not possible to apply the Jordan- Wigner transformation directly to (A.l) 
because of the non-local term Ri3{u) in the Yang-Baxter equation (A.2). For this reason, 
we consider an equivalent i?-matrix Rjk{u) 

Rjk{u) = PjkRjkiu) 

^a(«)i±|^ + c(tx)i^|^ + 6(«)(a+a, + a7a+) 

+ d{u){a+a+ + aja^), (A.4) 
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which satisfies the Yang-Baxter equation in the form 

Ru{u - v)R23{u)Ru{v) = R23{v)Ru{u)R23{u - v). (A.5) 
Here Pj^ is the permutation matrix 

Pjk = ^ + (T+a; + a- at. (A.6) 

Now we can apply the Jordan- Wigner transformation (A. 3) to (A. 4) and (A.5). Then 
we obtain a fermionic solution 

Ujk{u) ^a{u){njnk + (1 - nj){l - Uk)} + c{u){nj{l - Uk) + (1 - nj)nk} 

+b{u){c]ck + clcj) + d{u){c]cl - CjCk), (A.7) 

satisfying 

7^12(^^ - V)n23{u)ni2{v) = n23{v)ni2{u)n23{u - V). (A.8) 

In the derivation, the following transformation laws are of particular importance, 

0-7(77+1 — > -CjCj+i. (A.9) 

Strictly speaking, TZj^iu) is defined in (A.7) only for k = j + 1. However we extends its 
definition to arbitrary j and k by the formula (A.7). 

Finally multiplying the products of the fermionic permutation operators (2.7) 

^2^3^23 = ^23^13^12, (A. 10) 

to (A.7) from the left, we obtain 

n^2{u - v)ni3{u)n23{v) = 7^23(^^)7^13(«)7^12(M - v), (a.ii) 

with 

T^jk{u) = VjkTljkiu) 

= a{u){-njnk + (1 - nj){l - Uk)} + b{u){nj{l - n^) + (1 - nj)nk} 

+ c{u){c]ck + clcj) - d{u){c]cl + c,Ck). (A.12) 

This is identical to the fermionic i?-operator (2.2). In this way we can obtain the fermionic 
i?-operator (2.2) from the i?-matrix of the XY Z spin model. 

Just in the same way, the fermionic /^-operator for the ID Hubbard model can be 
obtained from the i?-matrix of the coupled spin model (1.5). 
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